Process reconstruction: Prom unphysical to physical maps via maximum likelihood 
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We show that the method of maximum likelihood (MML) provides us with an efficient scheme 
for reconstruction of quantum channels from incomplete measurement data. By construction this 
scheme always results in estimations of channels that are completely positive. Using this property 
we use the MML for a derivation of physical approximations of un-physical operations. In particular, 
we analyze the optimal approximation of the universal NOT gate as well as a physical approximation 
of a quantum nonlinear polarization rotation. 



Any quantum dynamics [1,2], i.e. the process that 
is described by a completely positive (CP) map of a 
quantum-mechanical system can be probed in two dif- 
ferent ways. Either we use a single entangled state of 
a bi-partite system [3] , or we use a collection of linearly 
independent single-particle test states [4,1] (that form 
a basis of the vector space of all hermitian operators). 
Given the fragility of entangled states in this Letter we 
will focus our attention on the process reconstruction us- 
ing only single-particle states. 

The task of a process reconstruction is to determine an 
unknown quantum channel (a "black box" ) using correla- 
tions between known input states and results of measure- 
ments performed on these states that have been trans- 
formed by the channel. The linearity of quantum dynam- 
ics implies that the channel £ is exhaustively described 
by its action Qj — > g'j — £[Qj\ on a set of basis states, i.e. 
a collection of linearly independent states Qj , that play a 
role of test states. Therefore, to perform a reconstruction 
of the channel £ we have to perform a complete state to- 
mography [1] of g'y The number of test states equals d? , 
where d = dim Ti. is the dimension of the Hilbert space 
associated with the system. Consequently, in order to re- 
construct a channel we have to determine d^{d^ — 1) real 
parameters, i.e. 12 numbers in the case of qubit (c? = 2). 

In what follows we will assume that test states can be 
prepared on demand perfectly. Nevertheless, the recon- 
struction of the channel £ can be affected by the lack 
of required information due to the following reasons: i) 
each test state is represented by a finite ensemble, corre- 
spondingly, measurements performed at the output can 
result in an approximate estimation of transformed test 
states; ii) the set of test states is not complete; and Hi) 
incomplete measurements on transformed test states are 
performed. In these cases some of the parameters that 
determine the map £ cannot be deduced perfectly from 
the measured data. In order to accomplish the channel 
reconstruction additional criteria have to be considered. 

In this Letter we will pay attention to the case i), that 
is typical for experiments - one cannot prepare an infinite 
ensemble, so frequencies of the measured outcomes are 
only approximations of probability distributions. Conse- 
quently, the reconstruction of output states g'^ can lead us 
to unphysical conclusions about the action of the quan- 



tum channel. As a result we can find a negative operator 
g'p or a channel £, which is not CP. 

In what follows we will introduce and compare two 
schemes how to perform a channel reconstruction with 
insufficient measurement data. Firstly, we will consider 
a rather straightforward "regularization" of the recon- 
structed unphysical map. Secondly we will exploit the 
method of maximum likelihood (MML) to perform an 
estimation of the channel. We will use these methods 
to perform a reconstruction of maps based on numerical 
simulation of the anti-unitary universal NOT operation 
(U-NOT) [5]. This is a finear, but not a CP map and 
we will show how our regularization methods will result 
in optimal physical approximations of the U-NOT oper- 
ation. We will conclude that the MML is a tool that pro- 
vide us with approximations of non-physical operations. 
In order to demonstrate the power of this approach we 
will also apply it to obtain an approximation of a non- 
linear quantum-mechanical map, the so called nonlinear 
polarization rotation (NPR) [6]. 

Structure of qubit channels. Quantum channels are de- 
scribed by linear trace-preserving CP maps £ defined on 
a set of density operators [1,2,7]. The complete positivity 
is guaranteed if the operator Q.£ = £ ® 'I\P+\ is a valid 
quantum state [P+ is a projection onto a maximally en- 
tangled state]. Any qubit channel £ can be imagined as 
an affine transformation of the three-dimensional Bloch 
vector r (representing a qubit state), i.e. r — > = Tr+t, 
where T is a real 3x3 matrix and t is a translation [7]. 
This form guarantees that the transformation £ is her- 
mitian and trace preserving. The CP condition defines 
(nontrivial) constraints on possible values of involved pa- 
rameters. In fact, the set of all CP trace-preserving maps 
forms a specific convex subset of all affine transforma- 
tions. 

The matrix T can be written in the so-called 
singular-value decomposition, i.e. T ~ Ri/DRy with 
Ru,Rv corresponding to orthogonal rotations and D = 
diag{Ai, A2, A3} being diagonal where Afc are the singu- 
lar values of T . This means that any map £ is a member 
of less-parametric family of maps of the "diagonal form" 
i.e. £[g\ = U^elVgV^W where U,V are unitary 
operators. The reduction of parameters is very helpful, 
and most of the properties (including complete positiv- 
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ity) of £ is reflected by the properties of $g. The map 
£ is CP only if <I>f is. Let us note that ^£ is determined 
not only by the matrix D, but also by a new translation 
vector r = Rjjt, i.e. under the action of the map $£ the 
Bloch sphere transforms as follows rj — > r'^ — XjVj + tj. 

A special class of CP maps are the unital maps, that 
transform the total mixture into itself. In this case 
t = f = 0, and the corresponding map <&£ is uniquely 
specified by just three real parameters. The positivity of 
the transformation results into conditions \Xk\ < 1- 
On the other hand, in order to fulfill the CP condition 
we need that the four inequalities |Ai ± A2I < |1 ± A3I 
are satisfied. These conditions specify a tetrahedron ly- 
ing inside a cube of all positive unital maps. In this 
case the extreme points represent four unitary transfor- 
mations I,ax,cry,az (see Fig.l). 
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FIG. 1. Unital CP maps are embeded in the set of all pos- 
itive unital maps (cube). The CP maps form a tetrahedron 
with four unitary transformations in its corners (extremal 
points) I,x,y,z corresponding to the Pauli (T-matrices. The 
un-physical U-NOT operation (Ai = A2 = A3 = — 1) and its 
optimal completely positive approximation quantum univer- 
sal NOT gate (Ai = A2 = A3 = — 1/3) are shown. 

Qubit channel estimation and its regularization. Re- 
constructions of states and processes share many fea- 
tures. Therefore we briefly remind us basic concept of 
the state reconstruction using finite ensembles of identi- 
cally prepared states. In this case one can obtain from 
estimated mean values of observable a negative density 
operator of a qubit g — i(/ + r ■ a) with |r| > 1. The 
reconstructed operator has always unit trace, but the as- 
sociated vector r can point out of the Bloch sphere. One 
can argue that the proper physical state is the closest one 
to the reconstructed operator, i.e. a pure state with 
pointing into the same direction. Formally it corresponds 
to a multiplication r by some constant fc, i.e. Tc = kr. 
The correction by k can be expressed as 

g^ = kg+{l~k)]^I = ]^{I + kr-a) (1) 

and it can be understood as a convex addition of the to- 
tal mixture ^/ represented by the center of the Bloch 
sphere, i.e. = (0,0,0). In other words, the correc- 
tion consists of the addition of completely random and 



equally distributed events (clicks) to outcome statistics, 
i.e. addition 0/ random noise. 

As we have seen from above, an important role in the 
state reconstruction is played by the total mixture i/, 
which is an average over all possible states. An average 
over all CP maps is the map A that transforms the whole 
state space into the total mixture, i.e. A[g\ — \l [8]. 

The reconstruction of qubit channels consists of the 
known preparation of (at least) four linearly independent 
test states gj and a state reconstruction of correspond- 
ing four output states g'y The process estimation based 
on the correlations gj — > g'^ is certainly trace-preserving 
(if all g'j are positive). Though the complete positivity 
is problematic. The average channel A can be used to 
correct ("regularize") improper estimations £ to obtain 
a CP qubit channel £c 

£.^k£Hl-k)A^[l/^^). (2) 

This method of channel regularization uses the same 
principle as the method for states, i.e. it is associated 
with an addition of random noise into data. 

Let us try to estimate what is the critical value of fc, 
i.e. the amount of noise that surely corrects any positive 
map. Trivially, it is enough to set fc = 0. In this case we 
completely ignore the measured data and the corrected 
map is A. However, we are interested in some nontrivial 
lower bound, i.e. in the largest possible value of k that 
guarantees the complete positivity. Let us consider, for 
simplicity, that the map under consideration is unital. 
Then the worst case of a positive map that is not CP, is 
represented by the universal NOT operation. 

Universal NOT. The logical NOT operation can be 
generalized into the quantum domain as a unitary trans- 
formation |0) |1), |1) — » |0). However, this map is ba- 
sis dependent and docs not transform all qubit states 
into their (unique) orthogonal complements Such 
universal NOT {£mi : IV') ^ IV'J-)) is associated with the 
inversion of the Bloch sphere, i.e. r — > — r, which is not 
a CP map. It represents an unphysical transformation 
specified by Ai = A2 A3 = —1. The distance (see Fig.l) 
between this map and the tetrahedron of completely pos- 
itive maps is extremal, i.e. it is most un-physical map 
among linear transformations and can be performed only 
approximatively. A quantum "machine" that optimally 
implements an approximation of the universal NOT has 
been introduced in Ref. [5]. The machine is represented 
by a map £mt — diag{l, — 1/3, — 1/3, — 1/3}. The dis- 
tance [8] between the U-NOT and its optimal physical 
approximation reads 

^(^NOT, ^not) = / dgTrK^NOT - ^not)!^]! = 1/3 . (3) 

J states 

The CP conditions imply that the minimal amount 
of noise necessary for a regularization of the universal 
NOT gate corresponds to the value of fc = 1/3, i.e. 
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Ai = A2 = As = —1/3 (see Fig.l). The channel repre- 
senting this point corresponds to the best CP approxima- 
tion of the universal NOT operation, i.e. to the optimal 
universal NOT machine originally introduced in Ref. [5]. 

One way how to interpret the "regularization" noise is 
to assume that the qubit channel is influenced by other 
quantum systems (the physics behind a dilation theorem 
[5]). 

The reason why we have to consider a noise in a recon- 
struction of quantum maps is that we deal with incom- 
plete measurement statistics (e.g., test states are rep- 
resented by finite ensembles). As a result, the recon- 



structed assignment gj 



— £[gj] is determined not 



only by the properties of the map £ itself but also by 
the character of the estimation procedure. In this situa- 
tion, the map itself can be unphysical, but if we request 
that the estimation procedure is such that the complete 
positivity of the estimated map is guaranteed then the 
result of the estimation is a physical approximation of an 
unphysical operation. In order to proceed we assume the 
method of maximum likelihood. 

The maximum likelihood method is a general estima- 
tion scheme [9,10] that has already been considered for 
a reconstruction of quantum operations from incomplete 
data. It has been studied by Hradil and Fiurasek [11], 
and by Sachci [12] (criticized in Ref. [13]). The task of 
the maximum likelihood in the process reconstruction is 
to find out a map £, for which the likelihood is maximal. 
By the definition we assume that the estimated map has 
to be CP. Let us now briefly describe the principal idea 
in more details. 

Given the measured data represented by couples gk , Fk 
{gk is one of the test states and Fk is a positive operator 
corresponding to the outcome of the measurement used in 
the kth run of the experiment) the likelihood functional 
is defined by the formula 
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L{£) = -\ogY[pik\k) = -J2^ogTT£[gk]Fk 



(4) 
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where N is the total number of "clicks" and we used 
p{j\k) — TT£[gk]Fj for a conditional probability of using 
test state gk and observe the outcome Fj. The aim is 
to find a physical map £est that maximizes this function, 
i.e. L(£esf ) = maxg L(£'). This variational task is usually 
performed numerically. 
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FIG. 2. The distance d{£mT, Sest) as a function of the num- 
ber of measured outcomes A'^ in the logarithmic scale. We 
used 6 input states (eigenvectors of axjO'yjaz) and we "mea- 
sured" operators ax, cry, a^- The distance converges to the 
theoretical value 1/3 that corresponds to the optimal univer- 
sal NOT. 



0.! 
0.6 
0.4 
0.2 ■ 


-0.2 







'^(a) ^^^^^ 




f . , 


} I 





10 12 14 16 



FIG. 3. We present analytical as well as numerical results 
of an approximation of a non-linear map £g given by Eq. (6) 
for different values of the parameter 6. The numerical ("ex- 
perimental") results shown in the graph in terms of a set of 
discrete points with error bars are obtained via the MML. For 
every point (9), the nonlinear operation was applied to 1800 
input states that have been chosen randomly (via a Monte 
Carlo method). These input states have been transformed 
according to the nonhnear transformation (6). Subsequently 
simulations of random projective measurements have been 
performed. With these "experimental" data a linear opera- 
tion was numerically searched for, which maximizes the prob- 
ability to obtain the same results. The resulting approxima- 
tion specified by a value of A (error bars shown in the graph 
represent the variance in outcomes for subsequent runs) trans- 
forms the original Bloch sphere as it is shown in an inset for 
a particular value 9 — 3. The figure (a) corresponds to result 
obtained by MML, and the figure (b) has been obtained via 
analytic calculations. We see that the original Bloch sphere is 
transformed into an ellipsoid, one axis of which is significantly 
longer than the remaining two axes, that are of a compara- 
ble length. The mean of these two lengths corresponds to 
the parameter A that specifies the map. The analytical ap- 
proximation £e of the nonlinear NPR map is characterized by 
the parameter A that is plotted (solid line) in the figure as a 
function of the parameter 9. 



Numerical results. Our approach is different from 
those described in Rcfs. [11-13] in the way, how we find 
the maximum of the functional defined in Eq.(4). The 
parametrization of £ itself, as defined by Eq.(2), guar- 
antees the trace-preserving condition. The CP condition 
is introduced as an external boundary for a Nelder-Mead 
simplex scheme. Therefore, there is no need to intro- 
duce Lagrange multipliers. We use the eigenstates of 
ax,cry,az as the collection of test states. The data are 
generated as (random) results of three projective mea- 
surements axjCTy, az applied in order to perform the out- 
put state reconstruction. In order to analyze the con- 
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vergcncc of the method we have performed the recon- 
struction for different number of chcks and compare the 
distance between the original map £hot and the estimated 
map £est- The result is plotted in Fig. 2, where we can 
see that the distance converges to 1/3 as calculated in 
Eq.(3). For N = 100 x 18 clicks, i.e. each measurement 
is performed 100 times per particular input state, the 
algorithm leads us to the map 



1 
-0.0002 -0.3316 
0.0138 -0.0031 
-0.0137 0.0298 




-0.0074 0.0203 
-0.3334 0.0488 
-0.0117 -0.3336 



(5) 



which is very close [d{£est,£app) = 0.0065] to the best 
approximation of the NOT operation, i.e. £app = 
diag{l,-l/3,-l/3,-l/3}. 

We conclude that for large N the MML reconstruction 
gives us the same result as a theoretical prediction de- 
rived in Ref. [5]. ^From here it follows that the MML 
helps us not only to estimate the map when just in- 
complete data are available, but also serves as a tool 
to derive physical approximations of unphysical maps. 
The reconstruction procedure guarantees that the esti- 
mation/approximation is physical. In order to illustrate 
the power of this approach we will find an approxima- 
tion of a non-linear quantum mechanical transformation 
that is even "more" unphysical than the linear though 
antiunitary U-NOT operation. 

Nonlinear polarizatAon rotation Let us consider a non- 
linear transformation of a qubit defined by the relation 
[6]: 



ge 



(6) 



Unlike the universal NOT this map is nonlinear. Four 
test states are not sufficient to get the whole information 
about the action of nonlinear maps. Consequently, the 
fabricated data must use all possible input states (that 
cover the whole Bloch sphere) as test states, but still we 
use only three different measurements performed on out- 
comes that are sufficient for the state reconstruction. We 
note that a straightforward regularization via the addi- 
tion of noise cannot result in a CP map unless the orig- 
inal map is not completely suppressed by the noise, i.e. 
the regularization leads to a trivial result £ = A. How- 
ever, as we shall see, the maximum likelihood approach 
gives us a reasonable and nontrivial approximation of the 
transformation (6). 

Firstly, we present an analytic derivation of a physical 
approximation of £e. This approximation is the closest 
physical map £g, i.e. d{£g,£o) = min. The map £g ex- 
hibits two symmetries: the continuous U(l) symmetry 
(rotations around the ^-axis) and the discrete ax sym- 
metry (rotation around the a;-axis by tt). The physical 
approximation £9 should possess these properties as well. 
Exploiting these symmetries the possible transformations 
of the Bloch vector are restricted as follows x — » Xx, y 



Ay, z pz. In the process of minimalization the param- 
eter p behaves trivially and equals to one. It means that 
£g is of the form £\ = diag{l. A, A, 1}. Our task is to 
minimize the distance d{£g,£\) = J dg\£g[g] — £\[g] \ and 
to find the physical approximation £g, i.e. the functional 
dependence of A on 6. 

We plot the parameter A that specifies the best physi- 
cal approximation of the NPR map in Fig. 3. In the same 
figure we also present a result of the maximum likelihood 
estimation of the NPR map based on a finite number of 
"measurements" (for technical details see the figure cap- 
tions). We conclude that the MML is in an excellent 
agreement with our analytical calculations. 

In summary, we have shown that theoretical methods 
that are designed for a reconstruction of quantum maps 
from incomplete date can be modified (see the method of 
maximum likelihood with an explicitly incorporated con- 
dition of the complete positivity) for an efficient deriva- 
tion of optimal (given the figure of merit) approxima- 
tions of un-physical maps. The derivation of physical 
approximations of un-physical maps based on reconstruc- 
tion schemes as discussed in the Letter is very natural. 
This view is siipported by the fact that numerical "recon- 
structions" are compatible with the analytical results. It 
should be pointed out that we have chosen two relatively 
simple examples (the universal NOT gate that is a non- 
CP linear map and the nonlinear polarization rotation) 
where analytical solutions can be found due to intrinsic 
symmetries of the considered maps. As soon as the maps 
do not posses such symmetries analytical calculations are 
essentially impossible and the power of our method be- 
comes obvious. 
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